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Abstract 

Liouville integrable separable systems with quadratic in momenta 
first integrals are considered. Particular attention is paid to the sys- 
tems generated by the so-called special conformal Killing tensors, i.e. 
Benenti systems. Then, infinitely many new classes of separable sys- 
tems are constructed by appropriate deformations of Benenti class sys- 
tems. 

1 Introduction 

The separation of variables for solving by quadratures the Hamilton-Jacobi 
(HJ) equations of related Liouville integrable dynamic systems with quadratic 
in momenta first integrals has a long history as a part of analytical me- 
chanics. There are some mile stones of that theory. First, in 1891 Stackel 
initiated a program of classification of separable systems presenting condi- 
tions for separability of the HJ equations in orthogonal coordinates ^-IS]- 
Then, in 1904 Levi-Civita found a test for the separability of a Hamilto- 
nian dynamics in a given system of canonical coordinates [3]. The next 
was Eisenhart 0-0) who in 1934 inserted a separability theory in the con- 
text of Riemannian geometry, making it coordinate free and introducing the 
crucial objects of the theory, i.e. Killing tensors. This approach was then 
developed by Woodhouse |H], Klanins and others. Finally, in 1992, 

Benenti - JHj constructed a particular but very important subclass of 
separable systems, based on the so called special conformal Killing tensors. 

The first constructive theory of separated coordinates for dynamical sys- 
tems was made by Sklyanin ^Ij. He adopted the method of soliton theory. 
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i.e. the Lax representation, to systematic derivation of separated coordi- 
nates. In that approach involutive functions appear as coefficients of char- 
acteristic equation (spectral curve) of Lax matrix. The method was success- 
fully aplied to separation of variables for many integrable systems |14|-jl9|. 

Recently, a modern geometric theory of separability on bi-Poisson mani- 
folds was constructed [in]-HHl) related to the so-called Gel'fand-Zakharevich 
(GZ) bi-Hamiltonian systems |2ni)|23|- Obviously, it contains as a spe- 
cial case Liouville integrable systems with all constants of motion being 
quadratic in momenta functions. Indeed, Ibort et.al. |i2^ proved that the 
Benenti class of systems can be lifted to the GZ bi-Hamiltonian form. 

In the following paper we construct in a systematic way all Liouville in- 
tegrable systems on Riemannian spaces, which are of the GZ type, including 
as a special case the Benenti class of systems. What is important, infinitely 
many classes of separable systems are constructed from appropriate defor- 
mations of the Benenti class of systems. In that sense we demonstrate the 
crucial role of this particular class of systems in the separability theory of 
dynamic systems with quadratic in momenta first integrals. 

2 Separable dynamics on Riemannian spaces 

Let {Q,g) be a Riemann (pseudo-Riemann) manifold with covariant met- 
ric tensor g and local coordinates q^,...,q^. Moreover, let G := be a 
contravariant metric tensor satisfying Yl^=i9ij^''^ ~ ^i- "^^^ equations 

qit + ^kQiQt=G"'dkV{q), i=l,...,n, ^ ^ (1) 

describe the motion of a particle in the Riemannian space with the metric 
g, where L*-^ are the Levi-Civita connection components. Eqs.(^ can be 
obtained by varying the Lagrangian 

C{q,qt) = ^^g^jqlqi-V{q) (2) 

hi 

and are called Euler-Lagrange equations. Obviously, for G = / equations 
reduce to Newton equations of motion. 

One can pass in a standard way to the Hamiltonian description of dy- 
namics, where the Hamiltonian function takes the form 

i=l * «J=1 j 
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and equations of motion are 




i _dH _ dH 

Pj, — \-I J \§g J ''''~~W 

Xh denotes the Hamiltonian vector field with respect to a canonical Poisson 
tensor ^tnd the whole dynamics takes place on the phase space J\f = T*G 
in local coordinates {q^, ...,Pn)- 

Of special importance is the geodesic motion V{q) = 0, with Euler- 
Lagrange equations and Hamiltonian representation in the form 

qit+rikQkt =0, i = l,...,n 

t 

n \ 1 " 

^ j =eodE = XE, E=-Y,G^^PiPr 

Functionally independent Hamiltonian functions Hi^i = 1, ...,n are said 
to be separable in the canonical coordinates (A, (£) if there are n relations, 
called the separation conditions (Sklyanin (14j). of the form 



(pi{X\Hi;Hi,...,Hn) = 0, i = l,...,n, det 



d(pi 



+ 0, (3) 



which guarantee the solvability of the appropriate Hamilton-Jacobi equa- 
tions and involutivity oi Hi. A special case, when all separation relations 
(jSJ are affine in ffj, is given by the set of equations 

n 

^Y^^\[\i,[ii)B.k = '4)i{\i,\Xi), i = l,...,n, (4) 

fc=i 

where and are arbitrary smooth functions of their arguments, is called 
the Stdckel separation conditions and the related dynamic systems are called 
Stackel separable. 

We are going to present a subclass of one-particle dynamics, containing 
Liouville integrable and separable systems with n quadratic in momenta 
constants of motion. Stackel ^-[21 was a first who gave the characterization 
of equations of motion integrable by separation of variables. He proved that 
if in a system of orthogonal coordinates (A, /i) there exists a non-singular 
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matrix 93 = ('/'^(Afc)), called a Stdckel matrix such that the Hamiltonians 
Hj. are of the form 

Hr = \j2^^-'yr{^^^+a.{X^)), (5) 
i=l 

then Hr are functionally independent, pairwise commute with respect to the 
canonical Poisson bracket and the Hamilton-Jacobi equation associated to 
Hi is separable. Indeed, for quadratic in momenta constants of motion, the 
Stackel separation conditions @ take the general form 

n 

Y^c|>H^')Hk = -f^{X')^^^+7^{X') z = l,...,n, (6) 

k=l 

where fi,ji, (f)^ are arbitrary smooth functions of its argument and the nor- 
malization (j)^ = 1, i = l,...,n is assumed. To get the explicit form of 
Hk = i/fc(A,^) one has to solve the system of linear equations ©• The 
results are the following 



' /i(Ai) /i(Ai) /i(Ai) /i(Ai) \ 



\ <^i(A") 0?,(A") 0r'(A") 1 ; 



Eisenhart considered Stackel separable systems in the frame of one- 
particle dynamics on Riemannian (pseudo-Riemannian) space {Q,g)- He 
gave a coordinate-free representation for Stackel geodesic motion introduc- 
ing a special family of Killing tensors [S]-[7]. As known, a (1, l)-type tensor 
K = (Kj) (or a (2, 0)-type tensor KG = A = {A'^)) is called a Kilhng tensor 
with respect to g if 



where {■, ■}eo means a canonical Poisson bracket. He proved [3]- [7] that the 
geodesic Hamiltonians can be transformed into the Stackel form @ if the 
contravariant metric tensor G = g~^ has (n — 1) commuting independent 
contravariant Killing tensors Ar of a second order such that 



Er = ^^Ai^piPj, 
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admitting a common system of closed eigenforms Oi 

{A* — v],G)ai = 0, dai = 0, i = 1, n, 

where v). are eigenvalues of (1, 1) Killing tensor = A^g {K* = gA*). 

For n degrees of freedom, let us consider n Stackel Hamiltonian functions 
in separated coordinates in the following form 

n 
i=l 

where // = (/ii, ^„)^and Vr{\) are appropriate potentials separable in 
(A, ^) coordinates. For the integrable system Q calculated from © 

,^,fj{}^)d^W^ +idetW^ 



,..det W^^'- 

i=l 



where 



1=1 



w 




and W^^ is the (n — 1) x (n — 1) matrix obtained from W after we cancel its 
ith row and fcth column. 

In our further considerations we restrict to the so-called GZ case, when 
(^f(Ai) are monomials, fi = f and ji = 7. Then, the Stackel separation 
conditions are n copies of a so-called separation curve 

H^C + - + HnC" = ^/(e)/^' + 7(e), m„ = < m„_i < ... < mi G N, 

(8) 

with = (A%/ii), i = 1, ...,n. 

In (A, /u) coordinates n related Hamilton-Jacobi (HJ) equations 

dW 

for a generation function W{X,a) = ^j(Aj, a), decouples into n ordi- 

nary differential equations 

1 /dW,^ ^ 
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and hence, the imphcit solution of dynamical system with Hamiltonian func- 
tion Hr is given by 



E 

k=l 



= trSri + consti, i = n, 



where 2/(^)[a(^) — 7(1^)] = called the inverse Jacobi problem. 

In this context, a question about classification and construction in nat- 
ural coordinates of all separable systems on Riemannian spaces, with n 
quadratic in momenta constants of motion, arises. The classification can be 
made with respect to the admissible forms of Stackel separability conditions. 
The right hand side of the conditions © is always the same for the class of 
systems considered 

r.h.s. = ^ ^ i;{X\i^,), (9) 

so different classes of separable systems are described by different forms 
of the l.h.s. of Stackel conditions, while systems from the same class are 
described by different / and 7 in @. 



3 Separable systems in natural coordinates 

Among all Stackel systems a particularly important subclass consists of these 
considered by Benenti , jT^l , ^3] and constructed with the help of the 
so-called special conformal Killing tensor. Let L = (L* ) be a second order 
mixed type tensor on Q and let L : Af — > M be a function on M defined as 

{L,E}g^^=KE, where K = {e,E}gg, e = Tr{L), 

then L is called a special conformal Killing tensor with the associated poten- 
tial £ = Tr{L). An important property of L is the vanishing of its Nijenhuis 
torsion. 

For the Riemannian space (Q,g,L), geodesic flow has n constants of 
motion of the form 

Er = -Y. Al^p^pj = - ^ {KrGfpip,, r = 1, n, (10) 

i,j=l ij=l 
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where and Kr are Killing tensors of type (2,0) and (1,1), respectively. 
Moreover, as was shown by Benenti ^^,^21) all the Killing tensors Kr with 
a common set of eigenvectors, are constructed from L by 

r 
k=0 

where pr are coefficients of the characteristic polynomial of L 

deti^I - L) = C + PiC'^ + - + Pn, Po = l, (11) 
or equivalently by the following 'cofactor' formula j29|.[3U| 

n-l 

C0f{(I-L) = Y,Kn-^e, 
1=0 

where cof{A) stands for the matrix of cofactors, so that cof{A)A = (det A)I. 
So, for a given metric tensor g, the existence of a special conformal Killing 
tensor L is a sufficient condition for the geodesic fiow on Q to be a Liouville 
integrable Hamiltonian system with all constants of motion quadratic in 
momenta. Moreover, the basic separable potentials Vr"^^ are given by the 
following recursion relation [201, 1221 

y^(-+i)=p,y/-)-y^(-), (12) 

where the first nontrivial potentials are 

It turns out that with the tensor L we can (generically) associate a co- 
ordinate system on Q in which the geodesic flows associated with all the 
functions Er separate. Namely, let {\^{q), X"{q)) be n distinct, function- 
ally independent eigenvalues of L, i.e. solutions of the characteristic equation 
det(^/ — L) = 0. Solving these relations with respect to q we get the trans- 
formation X ^ q : qi = ai{\). The remaining part of the transformation to 
the separation coordinates can be reconstructed from the generating func- 
tion W{p,\) = ^iPiOii{X). In the {X, p) coordinates the Stackel separation 
conditions ^ for Benenti Hamiltonian functions Hr = Er + Vr are 
given by the separation curve of the form 

HiC-' + H^C-^ + ... + Hn = \f{i)p^ + j = 0, 1, 2, ... . (13) 

It is important to notice, that all Liouville integrable systems of classical 
mechanics, with quadratic in momenta first integrals, that was separated in 
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XIX and XX centuries, belong to the Benenti class. The reason is that only 
the Benenti class contains dynamic systems on flat and constant curvature 
Riemannian spaces. Another important feature of Benenti systems is that 
all of them can be lifted to a one Casimir bi-Hamiltonian form (28j|. 

Now we present how to construct all remaining classes of separable sys- 
tems by an appropriate deformations of the Benenti class. Let start from 
the separability condition Q for n Hamiltonian functions in the following 
form 



1 



n+j 



where H. 



ni 



"2 



/c G N, j >k, 
(14) 

0, 1 < ni < ... < n^, < n + /c, and the 



separability condition for Benenti systems with the same ip 



(15) 



Theorem 1 Deformation of the Benenti system with separation curve \15\) 
to the new system with separation curve is given by a following deter- 
minant form 



Hr-k Pr-1 
Hfii—k Pni — l 

^nu—k Pn^—l 



Pr-k 
Pn\~k 

Prik-k 



"1^ 



(16) 



where pi, i = 0, ...,n are coefficients of the characteristic polynomial of the 
conformal Killing tensor L related to the Benenti system. 

The formula (jlGj) applies separately to the geodesic and the potential 
parts. The geodesic part Er can be presented in the following form 



Er 



1 



p KrGp, r = l,...,n + k, 



(17) 



where metric tensor G and Killing tensors Kr are 



1 



G = {-if-DoG, Kr = Kr- Kr-lDiD^' + ... + {-lfKr-kDkDo\ 
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Pni~k 



Prik-l 



ni — 1 



Kni-k 



K. 



rii 



K.. 



n-i—i+l 



K. 



"1- 



ni — fc+1 



K 



rik-i+l 



K.. 



Tife— i— 1 



and Kfn in determinant calculations are treated as symbols not matrices. 
The proof of the theorem as well as other details the reader can find in |3f)j . 

Such constructed systems, although obtained through the deformation 
procedure on the level of Hamiltonian functions, are far from being trivial 
generalizations of Benenti systems. There is no obvious relations between 
solutions of a given Benenti system and all its deformations. In each case we 
have a different inverse Jacobi problem to solve. Notice, that the common 
feature of appropriate deformed systems is the same set of separated coordi- 
nates, determined by the related Benenti system. Moreover all of them can 
be lifted to multi-Casimir bi-Hamiltonian form |3Uj . 



4 Example 

Consider the case n = 2. Let Q be a two dimensional flat space parametrized 
by canonical coordinates q = {q^,q^) with the contravariant metric tensor 
G and related special conformal Killing tensor L of the form 




Two geodesic Hamiltonians Ei and E2, according to UlUj) are 

^1 = ^Pi + ^P2^ E2 = ^gVp2 - \'i^Pl- 
Let us choose the potential V^'^\ As pi = —q^ and p2 = —5(9^)^ hence 

V?^ = {q'f + \q\q')\ = ^^{q^ + lilW?- 

It is one of the integrable cases of the Henon-Heiles system with Hamiltonian 

(3) (3) 

function Hi = Ei+ , second constant of motion H2 = E2+ V2 and 
Newton equations 

{q% = -3{q'f - \{qr 
{q')u = -qV. 
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The transformation to separated coordinates (A, /x) takes the form 

AVi , A2/7,2 nTT2 f ^'1 , ^2 \ 

^^ = A^+A^' ^^ = ^-^^ Va^ + A^J' 
and the separation curve is 

Now, let us consider the simplest deformation given by A; = 1, ni = 2. 
Then, 

Vi = --V^i, ^^2 = 0, % = V2--Vu 

Pi Pi 

hence 

i^i and arc in involution and are separated in the same coordinates 
(A, ^) as the Henon-Heiles system. The appropriate separation curve takes 
the form 
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